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The WILLI detector, built in IFIN-HH Bucharest, in collaboration with KIT Karlsruhe, is a rotatable modular 
detector for measuring charge ratio for cosmic muons with energy < 1 GeV. It is under construction a mini-array 
for measuring the muon charge ratio in Extensive Air Showers. The EAS simulations have been performed with 
CORSIKA code. The values of the muon flux, calculated with semi-analytical formula, and simulated with 
CORSIKA code, based on DPMJET and QGSJET models for the hadronic interactions, are compared with 
the experimental data determined with WILLI detector. No significant differences between the two models and 
experimental data are observed. The measurements of the muon charge ratio for different angles-of-incidence, 
(performed with WILLI detector) shows an asymmetry due to the influence of magnetic field on muons tra- 
jectory; the values are in agreement with the simulations based on DPMJET hadronic interaction model. The 
simulations of muon charge ratio in EAS performed with CORSIKA code based on three hadronic interaction 
models (QGSJET2, EPOS and SYBILL) show relative small difference between models for H and for the Fe 
showers; the effect is more pronounced at higher inclination of WILLI detector. The future measurements should 
indicate which model is suitable. 



1. WILLI detector 

The rotatable WILLI detector [1], shown in FigJTJ 
is a modular system, each module being formed by 
a scintillator layer, 3 cm thickness in Al frame, 1 cm 
thickness, operating for all azimuth angles and down 
to 45° zenith inclination. 





Figure 1: Left: The vertical WILLI detector Right: The 
rotatable configuration of WILLI 

The detector permits the measurement of the en- 
ergy deposited by traversing muons and additional, 
for muons stopped in the detector, the moment of de- 
cay up to 50 [is. The active layers of the detector 
are plates of scintillators of 1 m 2 which are read by 
two photomultipliers mounted on opposite corners of 
each plate, using alternatively corners. The main aim 
of WILLI detector is to determine the muon charge 
ratio by measuring the life time of stopped muons in 
the detector layers: the stopped positive muons de- 
cay with a lifetime of 2.2 /is, while negative muons 
are captured in the atomic orbits, leading to an ef- 



fectively smaller lifetime depending on the stopping 
material. The muon charge ratio is determined from 
the measured decay curve of all muons stopped in the 
detector, by fitting the measured decay spectrum with 
the theoretical curve. 



2. The simulations of the cosmic muon 
flux 



The simulations of the EAS have been performed 
with CORSIKA program [2[, based on six differ- 
ent models for the description of the high - energy 
hadronic interaction: DPMJET II.5 [3], QGSJET [|, 
VENUS @, SIBYLL @ and three different models for 
the description of low - energy hadronic interaction: 
GHEISHA [7], UrQMD 1.1 [8] and DPMJET, which 
includes some extensions allowing the simulation of 
hadronic interaction down to 1 GeV; the threshold is 
set by default to Elab = 80 GeV/n. The local geomag- 
netic field is included in CORSIKA in the approxima- 
tion of homogeneous field, as described by the Interna- 
tional Geomagnetic Reference Field for the year 2000 
Q . The geomagnetic cutoff is calculated with a Monte 
Carlo procedure of the possible particle trajectories in 
the back-tracking method, enabling the calculation of 
a table of allowed and forbidden trajectories 1 1 Oil . The 
particle tracking is based on GEANT 3.21, [ll| start- 
ing at 112.83 km, the top of the atmosphere as defined 
in CORSIKA. 
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3. The investigation of the cosmic muon 
flux 

In view of detailing geomagnetic effects, the muon 
fluxes have been calculated [12] for two different loca- 
tions with different magnetic cutoff: Hiroshima (34° 
N, 132° E) with the geomagnetic cutoff 11.6 GV and 
Bucharest (44° N, 26° E) with geomagnetic cutoff of 
5.6 GV [10|. Fig |2] (left) compares for Bucharest and 
Hiroshima the muon flux calculated with CORSIKA 
with the semi-analytical formulae of Nash [13| and 
Geiser Q. 
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Figure 2: The muon flux data compared with Monte 
Carlo simulations and semi-analytical formulae for 
Hirosima and Bucharest (left); The azimuthal variation 
of muon flux measured by WILLI compared with 
CORSIKA simulations(right) 



For the investigation of the azimuthal dependence 
of the charge ratio of atmospheric muons, a series 
of measurements [12| has been performed on four 
azimuth directions of incidence of the atmospheric 
muons: North, East, South, West, (N, E, S, W) for 
muons with inclined incidence, mean value at 35° and 
mean incident energy 0.5 GeV/c. The results are dis- 
played in Fig |3] (right). 
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Figure 3: Energy dependence of the muon charge ratio, 
measured separately for East and West directions, with 
comparing experimental data and simulations (left); The 
azimuthal variation of the muon charge ratio (right) 

The good agreement between the measurements 
and simulations data indicates that CORSIKA code 
describes well the azimuth variation i.e. the East- 
West effect observed by WILLI. The Okayama group 
reported less pronounced azimuth dependence [16j for 
muon incident energies above 1 GeV. 



We see that at low muon energies ( < 1 GeV) the 
muon flux is significantly influenced by the geomag- 
netic field and dependent on the particular observa- 
tion site, and from the direction of muon incidence. 
The muon flux for incident energies < 1 GeV was 
compared with CORSIKA simulations using DPM- 
JET and QGSJET2 models for hadronic interactions, 
see Fig|2] (right). 

The simulated data for atmospheric muons ob- 
tained with CORSIKA code, using both DPMJET 
and QGSJET interaction models, are in good agree- 
ment with measured data of WILLI detector for the 
whole azimuth range and a mean zenith angles of 35°. 



4. The investigation of the charge ratio 
for the cosmic muons 

The charge ratio R M (/i + //i _ ) of the atmospheric 
muons provides a sensitive test of the simulation of the 
fluxes as different charges have different path lengths 
from the production level to the detector, and conse- 
quently the decay probability for low-energy muons is 
modified. This influence leads to the latitude effects 
of the flux and to the so-called East West effect [lEj . 
which is more pronounced at lower muon energies and 
with larger observation angles, see Fig |3] (left). 



5. The detector WILLI-EAS 

We are building a small scintillator detector array 
nearby WILLI in order to trigger the muon detection 
of WILLI by events of small showers and to determine 
the average charge ratio within EAS [ItJ- The mini- 
array is formed by 12 detector stations, each being a 
scintillator plate, 3 cm thickness, divided in 4 parts 
(0.475 x 0.475 m2 - see FiggJ). 
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Figure 4: WILLI-EAS detector (left); The unit of the 
mini-array for detecting the shower events (right) 

It is expected that like for the East -West effect of 
atmospheric muons, the mean charge ratio is affected 
by the geomagnetic field, especially when - like it is the 
case with WILLI - low energy muons are registered. 
The simulations of muon charge ratio performed with 
CORSIKA based on QGSJET model performed for 
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H, see FigjSJ and Fe showers with primary energy 
10 15 eV, show azimuthal variation, more pronounced 
with increasing distances of WILLI detector, from the 
shower core jl8j . 



fill 


A R =(15-20 m) 
■ R =(45-50 m) 




* R =(145-150 m) 






. , , . i , . , . i , . , , i , . , , i 


, / 


. , , 1 , , . , 1 , , . , 1 



50 100 150 200 250 300 350 
Azimuth (°) 

Figure 5: Variation of the charge ratio R M of the mean 
muon density distribution of proton-induced EAS of 
inclined showers (9 — 45°) incident from the North with 
the primary energy of 10 15 eV at various distances R 
from the shower axis 

The simulations studies of the possible configura- 
tions of the mini-array have been performed with 
CORSIKA based on QGSJET model for a mini ar- 
ray with detector stations placed at different distances 
from WILLI detector and with different distances be- 
tween the modules. The quality of the shower core 
reconstruction was checked for all the cases. Figj6] 
shows the reconstructed shower cores for 250 H inci- 
dent showers placed in one position for 30° incident 
angles (left) and the quality of reconstruction for 750 
H incident showers uniform distributed on the arrays 
surface for 20, 30 and 45° zenithal incidences (right) 

0. 




x(m) R(m) 

Figure 6: The reconstructed shower cores for 250 H 
incident showers placed in the same position (0,0) for 30° 
(left) and the quality of the reconstruction for 750 H 
incident showers uniform distributed on the arrays 
surface for 20, 30 and 45° zenithal incidences (right), of 
the real configuration of the mini-array 

By simulating different incident H and Fe show- 
ers, with different energies and with different an- 
gles of incidence, the influence of the hadronic in- 
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teraction model on the EAS muon charge ratio has 
been investigated. Three different interaction mod- 
els, QGSJET2, EPOS and SYBILL, have been used. 
Fig 13 (left and center) compares the azimuth variation 
of muon charge ratio for the three models 

Using QGSJET2 and DPMJET interaction models, 
the simulations have been performed for the real con- 
ditions of the experiment, see FigJTl (right), observ- 
ing a pronounced asymmetry with increasing angle- 
of-incidence. 



6. Conclusions 

i) . The simulations of the muon flux for low en- 
ergy (< 1 GeV), using two interaction models (DPM- 
JET and QGSJET2) have been compared with results 
of the measurements performed with WILLI detector, 
noticing good agreement for both, with very small dif- 
ferences. 

ii) . The rotatable system WILLI allows measure- 
ments of the charge ratio of muons with different 
angles-of incidence, being possible to determine the 
East- West effect of the Earth's magnetic field. The 
measurements with rotatable WILLI, inclined at 35°, 
show a pronounced East- West effect, in good agree- 
ment with simulations data using CORSIKA code 
with DPMJET model, and with the East- West effect 
found in neutrino measurements (20j . 

hi). The simulations of H and Fe showers in- 
cident on WILLI-EAS detector, with different en- 
ergies and with different angles of incidence, have 
been compared for three different interaction models, 
QGSJET2, EPOS and SYBILL. The results of simu- 
lations show some difference between the three models 
and also between H and Fe showers. 

iv). The simulations of EAS muon charge ratio 
performed for the real condition of the experiment 
WILLI-EAS, using two interaction models, DPMJET 
and QGSJET2, indicate small differences between 
them, which become more pronounced for higher in- 
clination of WILLI. 
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Figure 7: The dependence of the charge ratio on the azimuth position of WILLI, at 15-20 m (left) and 45-50 m (center) 
from the shower axis, for proton and iron induced showers, coming from South at 30° zenithal inclination, using 3 
different hadronic interaction models (QGSJET2, EPOS and SYBILL); The variation with the azimuthal position of 
WILLI (inclined with 45°) of the EAS muon charge ratio for 5000 incident showers with 10 13 < Eprimary < 10 15 eV at 
45-55 m distance from the shower core simulated with 2 hadronic interaction models: QGSJET2 - full line and 
DPMJET - dotted line, (right) 
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